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We have used deprojected radial density and temperature 
profiles of a sample of 16 nearby CF clusters observed with 
XMM-Newton to test whether the effervescent heating model 
can satisfactorily explain the dynamics of CF clusters. For each 
cluster we derived the required extra heating as a function 
of cluster-centric distance for various values of the unknown 
parameters M (mass deposition rate) and fc (conduction effi- 
ciency). We fitted the extra heating curve using the AGN ef- 
fervescent heating function and derived the AGN parameters L 
(the time-averaged luminosity) and ro (the scale radius where 
the bubbles start rising in the ICM). While we do not find any 
solution with the effervescent heating model for only one ob- 
ject, we do show that AGN and conduction heating are not co- 
operating effectively for half of the objects in our sample. For 
^) most of the clusters we find that, when a comparison is pos- 
PsJ ' sible, the derived AGN scale radius ro and the observed AGN 
jet extension have the same order of magnitude. The AGN lu- 
d ' minosities required to balance radiative losses are substantially 
^! lowered if the fact that the AGN deposits energy within a finite 
' volume is taken into account. For the Virgo cluster, we find that 
the AGN power derived from the effervescent heating model is 
in good agreement with the observed jet power. 
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1. Introduction 

The radiative cooling time of the intracluster medium (ICM) 
in the cores of many galaxy clusters is short enough for the 
plasma to radiate an amount of energy equal to its thermal 
energy in less than a billion years. Therefore the gas should 
gradually condense, be replaced by the surrounding material, 
and ultimately form a cooling flow (hereafter CP). Early es- 
timates of the mass deposition rates from X-ray imaging (~ 
10^ - lO-'Moyr"') have indicated that the cores of CP clusters 
should contain large amounts of cold gas, but searches for these 
large amounts of condensed baryons have not been successful 
(see Pabian 1994 for a review of the subject of CPs before the 
launch of the XMM-Newton and Chandra satellites). 

XMM-Newton and Chandra observations have greatly im- 
proved our understanding of CP clusters. Most important, X- 
ray spectra failed to detect the emission lines that dominate the 
emission from gas below 2 keV (Peterson et al.'2001' Tamura 
et al. 2001 Kaastra et al. 2001, Peterson et al. 2003) , there- 
fore ruling out the standard CP model. In addition, as shown by 
Peterson et al. ("3003 ), the observed spectra cannot be explained 
without any substantial modifications to the general CP pro- 
cess. Recent observations have also shown that mass deposition 
rates are significantly smaller than those estimated previously 
(Bohringer et al. 2001 David et al. 2001 Peterson et al. 2003 1. 
Even though this implies a less severe discrepancy between CP 
mass deposition and star-formation rates (McNamara2004"l, no 
agreement between the amount of cold baryons predicted from 
X-ray observations and the one observed in other wavelengths 
has been found yet. Spatially-resolved spectroscopy shows that 
the temperature in CP clusters drops towards the center to ap- 
proximatively one third of the cluster mean temperature, indi- 
cating that the gas is prevented from cooling below these cutoff 
temperatures (Peterson et al. 2001 Tamura et al. 2001 , Kaastra 
et al. l2004l . 

This recent evidence shows that the dynamics of the ICM 
in CP clusters is not solely governed by cooling of the ICM, 
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and that some heating and/or non X-ray cooling mechanisms 
must be investigated. While non X-ray, rapid cooling mecha- 
nisms have been studied (see Peterson et al. 120031 and refer- 
ences therein), the most appealing mechanisms involve heating 
processes. 

Most CP clusters host a central active galactic nucleus 
(AGN) with strong radio activity (Burns7990' Ball et al. ll993> 
and, most important, recent observations show that these ra- 
dio sources are interacting with the ICM and often displace the 
hot gas, leaving cavities in their wakes (e.g., see Blanton l2004l 
Birzan et al. 2004 and references therein). Hence, it has been 
realized that gas heating by the outflows of the central AGN 
can be a vital process for the dynamics of CPs. This heating 
mechanism involves buoyant plasma bubbles inflated by the 
AGN, which then rise through the cluster atmosphere, expand, 
and ultimately heat the ICM (Churazov et al. 2002 Briiggen 
& Kaiser 2001 Ruszkowski & Begelman 2002 Briiggen & 
Kaiser |2'002> . Although there is no consensus regarding the 
ability of AGN heating to prevent the formation of CPs and 
the efficiency of the processes by which energy is delivered to 
the ICM, there is evidence suggesting that the class of mod- 
els in which the AGN energy input alone balances radiative 
losses is unable to quench the CP (Brighenti & Mathews 2002. 
Zakamska & Naravan l2003> . 

Despite the fact that the presence of magnetic fields in clus- 
ters implies that thermal conductivity is only a fraction of the 
Spitzer value, thermal conduction by electrons might play an 
important role in CPs. While some theoretical work has pointed 
out that thermal conduction must be highly suppressed (e.g., 
see Pabian n"994> . recent theoretical papers show that conduc- 
tivity can be as high as a substantial fraction of the Spitzer rate 
(Narayan & Medvedev 2001. Gruzinov 2002) and therefore 
it has been reconsidered as a possible heat source candidate. 
However, it has been shown that heat conduction alone fails 
to balance radiative losses at the center of the clusters (Voigt 
et al. l2002l Z akamsk a & Naravan l2003l Voigt & Pabian l2004l 
Kaastra et al. "2004) and that models with heat conduction as 
the only heating source are unstable (Soker i2003l) . 
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Other, less popular heating mechanisms are: fluid tur- 
bulence generated by substructure motion or cluster merg- 
ers (Fujita et al. 2004 1, heating from intra-cluster supernovae 
(Domainko et al. .2004 1. heating by hadronic cosmic rays 
(Colafrancesco et al. P004 1, and heating by magnetic fields 
(Soker & SarazinT590 Makishima et al. 2001 1. 

The inability of heat conduction to balance radiative losses 
and the failure of some AGN heating models have motivated 
the development of models in which these two mechanisms co- 
operate. Models with heat conduction and AGN heating acting 
together are very attractive due to the complementary nature of 
the two processes: thermal conduction is effective/ineffective 
in the outer/inner regions of the cluster, and AGN heating 
is effective/ineffective in the inner/outer part. Ruszkowski & 
Begelman (2002i (hereafter RB02) have recently proposed 
such a model. The ICM density and temperature evolved ac- 
cording to their model reach a final stable configuration in 
agreement with the general shape of observed density and tem- 
perature profiles in CF clusters. RB02's model (or the effer- 
vescent heating model, hereafter) is the only proposed theoret- 
ical model that can be effectively tested against observations. 
Unfortunately this has only been done for the Virgo cluster 
(Ghizzardi et al. 2004), and it is not clear whether the model 
can give a satisfactory explanation of the dynamics of CF clus- 
ters in general. 

In the present work we address the latter question using a 
sample of CF clusters observed with XMM-Newton. Our sam- 
ple, which consists of 16 objects, is large for this kind of anal- 
ysis and all the objects were analyzed using an appropriate and 
homogeneous procedure. 

The layout of this paper is as follows. In Sect.|2lwe present 
the sample and briefly summarize the procedure used to derive 
the basic quantities used in our analysis. In Sect. |3 we give a 
description of RB02's model and present the equations used in 
Sect. 13 where we fit the model to the data. In Sect.|5lwe sum- 
marize the results of our analysis and present the conclusions. 

We adopt Hq = 70 kms ' Mpc"', Q.^ = 0.3 and Oa = 0.7 
in all the calculations throughout this paper 

2. Data 

The sample consists of the clusters listed in Table^] These ob- 
jects are part of the sample analyzed in Kaastra et al. (2004 
hereafter K04), who derived deprojected radial temperature and 
density profiles for all the clusters in the sample, including a 
description of the sample and an extensive presentation of the 
data analysis. From the K04's sample, we excluded 4 clusters: 
Coma, A 754 and A 3266, because these are non-CF clusters, 
and A 1837 since temperature and density could be measured 
in only one bin within the cooUng radius. Spectral fitting was 
done over the full 0.2-10 keV range and, in general, temper- 
atures and gas densities were computed for the innermost 8 
shells (see K04, Table 3 for the boundaries between the shells), 
except when the data became too noisy in the outermost shell. 
The sample used in the present work therefore consists of 16 
clusters with CF for which gas density and temperature radial 
profiles are well measured. Since thermal conduction depends 
on temperature gradients and the AGN heating in the RB02 
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model involves a dependency on ICM pressure, we model gas 
temperature and density profiles by fitting the data with ana- 
lytic functions. We used more that one fitting function for both 
quantities in order to quantify the difference due to the mod- 
elling. In the following we describe these functions and the fit- 
ting to the data. 

2.1. Modelling of gas temperature and density 

The deprojected temperature profiles of our clusters exhibit a 
self-similar shape studied in Piffaretti et al. J2005> : the tem- 
perature declines from the maximum cluster temperature at a 
break radius r\,i- moving outwards and shows the characteris- 
tic temperature decline towards the X-ray emission peak. Since 
we are interested in the central cooling region and the cooling 
radius for a cooling time of 15 Gyr, /?cooi is smaller than rbr for 
all the clusters, the temperature profiles can be simply modelled 
by a function that is monotonically raising with radius. Hence, 
for each cluster we select temperature bins inside the radius 
^T.max - '"bi and fit them using the following expressions: 

r(r)^ro + r, y y (i) 

or 

r(r) = f„-fiexp(-^). (2) 

In order to reduce the number of parameters here, we set 
T{r = 0) equal to the temperature of the central bin for both fits 
and use ju = 2 in Eq.[0(Allen et al. 2001). For Perseus, Virgo, 
A 262, A 496, and A 3 1 12, the temperature in the observed ra- 
dial range is monotonicafly raising with radius, and therefore 
^T.max is just the distance of the outermost bin from the cluster 
center The best-fit parameters and i?T,max are given in Tabled 
Both temperature parametrizations are used in the computation 
of thermal conduction and in the modelling of the gas pressure. 
Our main results presented in Sect. 0] are achieved using the 
parametrization in Eq. [fl and we use the second parametriza- 
tion given in Eq.|2lonly in order to explore the effect of a differ- 
ent modelling on our main results. The changes introduced by 
using Eq. |2l are, as we will show in Sect. 14.41 quite small and 
do not change the results obtained by using Eq.^ 

We model the gas density by using a single jS-model given 

by: 

p(r)=po[l+{r/r,fy\ (3) 

The density profile is fitted within /?p,max, which is the radius at 
the center of the last radial bin, where a robust estimate of gas 
density and temperature is possible (see Piffaretti et al. 120051 
for the bin selection criterion). The best-fit parameters for the 
single j6-model and ^p^max are given in Table |5] 

An alternative parametrization of the gas density profile is 
the more complex double y6-model, which is a popular general- 
ization of the single y6-model used to model the central surface- 
brightness excess observed in CF clusters. Unfortunately in this 
case, the gas density is modelled using the sum of two single 
j8-models, so the number of free parameters is doubled: poj, 
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Table 1. The radius at the center of the last radial bin considered in the temperature fits and the best-fit parameters from fitting 
Eqs.[2and |2]to the data (1 cr eiTors in parenthesis). Here, yu = 2 in Eq.[2is used and, for both fitting functions, T{r - 0) is set 
equal to the temperature of the central bin. 
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y6,, and re ,-, with i — 1,2. As a consequence, while fitting the 
single yS-model to the density profiles gives statistically signif- 
icant values for the best-fit parameters, the large number of 
parameters adopted in the double j6-model, together with the 
small number of bins in which the gas density is measured, do 
not allow a significant determination of the best-fit parameters. 
Therefore, we decided to present our main results in Sect. |3 
using the single /3-model. Nonetheless, in order to constrain to 
which extent the double y6-modelling changes the results, we 
also fit the density profiles using a double j6-model with a re- 
duced number of fitting parameters, by setting fix - jSi and 
fixing one of the core radii equal to the core radius rj (the core 
radius fj is also used as a supplementary choice). The derived 
profiles are then used to model the gas pressure as done by us- 
ing the single j6-model results. A discussion of the influence of 
this different modelHng on our main results is given in Sect. l4.4l 
below. 

3. Effervescent heating 

Assuming spherical symmetry, the thermodynamic equations 
describing the ICM can be written in the form; 



vpr = const : 



M 
4^ 



r^v dv 



M,„, = — — — 
ar 



kTr 



H = 
where 



G dr GfiMp 



dlnT dlnp 



+ 



d\nr d\nr 



pvkT 



3 d\\\T 
2 dlnr 



dlnp 



dlnr 



(4) 

(5) 
(6) 

(7) 



and V is the gas flow velocity which is taken positive outwards 
and M is the mass flow rate, which is therefore negative for an 
inflow and assumed, throughout the paper, to be constant. Here, 



Table 2. The radius at the center of the last radial bin con- 
sidered in the gas density fits and the best-fit parameters from 
fitting Eq.l^to the data (1 cr errors in parenthesis). 



Cluster 




Po 








(kpc) 


(10-^"g/m^^) 


(kpc) 




NGC 533 


160 


1.85(0.77) 


7(2) 


0.46(0.01) 


Virgo 


35 


8.14(0.85) 


3(1) 


0.33(0.01) 


A 262 


199 


4.54(0.91) 


5(1) 


0.33(0.01) 


Sersic 159-3 


499 


4.36(0.06) 


35(1) 


0.59(0.01) 


MKW9 


238 


67.7(16.29) 


19(8) 


0.32(0.02) 


2A 0335+096 


309 


8.24(0.49) 


20(2) 


0.53(0.01) 


MKW3s 


402 


2.61(0.12) 


35(2) 


0.49(0.01) 


A 2052 


319 


4.02(0.36) 


19(2) 


0.45(0.02) 


A 4059 


412 


1.71(0.18) 


48(7) 


0.50(0.02) 


Hydra A 


673 


5.13(0.70) 


29(5) 


0.52(0.02) 


A 496 


193 


5.14(0.14) 


16(1) 


0.42(0.01) 


A3112 


645 


4.15(0.35) 


41(4) 


0.55(0.02) 


A 1795 


553 


4.26(0.18) 


42(3) 


0.52(0.01) 


A 399 


606 


0.72(0.08) 


76(15) 


0.38(0.03) 


Perseus 


229 


7.27(0.52) 


28(3) 


0.53(0.02) 


A 1835 


831 


5.24(0.23) 


80(5) 


0.62(0.02) 



M,ot is the total gravitational mass within the radius r and p, T 
are the gas density and temperature, respectively (yu = 0.61 is 
the mean molecular weight). Further, s - n^MT) is the plasma 
emissivity and e* is the energy due to the inflow/outflow of the 
gas. Then, e is computed from the deprojected electron density 
He and A(r), the cooling function for a plasma with an average 
metallicity Z = 0.5 Z© losing energy by bremsstrahlung and 
line emission. The quantity H in the energy equation (Eq.|6} is 
thus an extra heating term. In the framework of the effervescent 
heating model developed in RB02, the extra heating is provided 
by buoyant bubbles injected into the ICM by the central AGN 
(see this section below). 
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The total gravitational mass M,oi given in Eq.|5]differs from 
the mass estimated under the assumption of hydrostatic equi- 
librium, since it includes a velocity term. For the values of M 
considered in Sect.l^below, this term is negligible, so the mass 
analysis presented in Piffaretti et al. (2005 1 gives accurate val- 
ues for the mass profiles of the clusters in the present work. 

For all the clusters in our sample, the entropy of the gas in- 
creases monotonically moving outwards, almost proportional 
to the radius (Piffaretti et al. I2005> . Given this evidence, the 
convection term in the original effervescent heating model is 
neglected, since a declining entropy profile is essential for con- 
vection to operate (see the discussion in Ghizzardi et al. 2004 1. 
The term Scond is the heating due to thermal conduction, which 
is given by 



^cond 



1 d 

dr 



dJ\ 



(8) 



where k is the conductivity. For an ionized plasma, k is the 
Spitzer conductivity: 



1.84 X 10-5 (r)5/2 _i _ _i 

K - Ks erg cm s K 

In A 



(9) 



where In A ~ 40 is the Coulomb logarithm. 

Magnetic fields at the jiG level are known to be present 
in clusters of galaxies, while higher values, up to tens of yuG, 
have been measured at the center of CF clusters (see Govoni 
& Feretti IT004I for a recent review). In the presence of such 
fields thermal conduction is suppressed below the Spitzer rate 
given in Eq. |5]by a factor ~ 100 - 1000 (Binney & Cowie 
119811 Chandran & Cowley 1998 1. This condition has been es- 
sential for the development of multiphase CF models (e.g., 
Fabian[r594), but recent work has shown that the level of sup- 
pression of thermal conduction might not be as high as previ- 
ously thought. In particular, Narayan & Medvedev (2001 1 have 
shown that if the magnetic field is chaotic over a wide range of 
length scales, thermal conduction is enhanced to ~ 1/5 times 
the Spitzer value. Gruzinov ( 2002 ) pointed out that the effective 
heat conduction in a random variable magnetic field is boosted 
to a factor of 3 below the Spitzer value. The latter result yields 
the upper limit for the efficiency of heat conduction in the ICM. 
Therefore, in our analysis presented in Sect. |4] we investigate 
models with thermal conduction varying from zero to the max- 
imum level of 1/3 times the Spitzer value. The ICM conductiv- 
ity K is therefore characterized by the fraction according to 

Since thermal conduction is most efficient at high temper- 
atures and it generally fails to balance radiative losses in the 
central parts of CF clusters, a heat source able to supply energy 
in the central parts is needed. Such a central heat source is in- 
corporated in the RB02 model, in which AGNs are assumed to 
inject buoyant bubbles into the ICM, which heat the ambient 
medium by doing pdV work as they rise and expand adiabat- 
ically. Here we summarize the derivation of the AGN heating 
function H'^'^'^ (Begelman i2001.) . Assuming a steady state, the 
energy flux available for heating is: 



where pt, is the pressure of the bubbles and yt, their adiabatic 
index. Assuming that the bubbles' partial pressure scales as the 
ICM pressure p, the heating function can be expressed 

according to: 



with 
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Idlnp 
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r d\nr 



[l-e-''''''')dr. 



(11) 



(12) 



(13) 



(10) 



Here, po is some reference ICM pressure (in the following we 
take its value at the cluster center) and L the time-averaged 
luminosity of the central source. The term 1 -exp(-r/ro) intro- 
duces an inner cutoff that fixes the scale radius where the bub- 
bles start rising in the ICM. Therefore an inner cutoff is already 
taken into account when performing the integration in Eq. [O] 
and with zero as the lower integration limit gives a correct es- 
timate for the time-averaged luminosity of the AGN. It must 
be noted, however, that the upper integration limit in Eq. [O] 
should be replaced by a finite number, i.e. the radius within 
which the bubbles effectively deposit heat into the ICM. This 
issue is explored in Sect. l4.5l since it is through Eq.[0]that the 
time-averaged luminosity of the central source is estimated and 
a finite integration limit has the effect of diminishing its energy 
requirement. 

Finally, we notice that the following results are not only 
valid for the AGN-heating mechanism just described (i.e., heat- 
ing by buoyant bubbles), but for any mechanism for which its 
dependence on the ICM pressure is equivalent to the one given 
in Eq.HI] 

4. Results 

From the observed deprojected gas temperature and density 
profiles, the required extra heating H can be computed for fixed 
values of the mass deposition rate M and the conduction ef- 
ficiency fc, and then fitted using the AGN heating function 
ffACN ^ For each bin with measured gas density and tempera- 
ture, we compute the gas emissivity s. For a fixed pair (M, fc), 
the conductive heat Scnd and the heat term s* are computed 
using the temperature and gas parametrizations given in Eqs.^ 
and since they depend on the gradients of these quanti- 
ties. Then Samd and e* are evaluated at the radius where the 
gas emissivity is computed, to finally obtain extra heating H 
through Eq.|6| 

The errors of H are computed from the errors (1 cr) of 
the gas emissivity and the ecnd and s* errors. The latter have 
been evaluated by propagating the 1 cr eiTors of the gas density 
and temperature best-fit parameters. While two different radial 
ranges have been used to fit the gas density and temperature 
profiles (^p,max and ^x.max, see Sect. 12. 1> . the extra heating H is 
computed in bins within /?T,max, since Rj.max < Rp.m-dx- Finally 
the AGN heating function given in Eq. [^is fitted to the extra 
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heating data points using a x~ minimization. The gas pressure 
and pressure gradients in Eq.f^are evaluated using Eqs.^and 
13 For any fixed pair (M, f^) the free parameters are therefore 
ro and the luminosity of the central AGN L. The bubbles' adi- 
abatic index yt, is fixed to 4/3 throughout the paper (i.e., we 
assume relativistic bubbles), since the inclusion of -yi, as an ad- 
ditional free parameter leads to very large errors in the derived 
best-fit parameters. A discussion of the effects due to the use of 
different gas density and temperature parametrizations is given 
in Section 

4.1. Zero AGN heating 

Despite the stability issues mentioned above, the case of ther- 
mal conduction by electrons acting as the only heat source has 
been already explored quite extensively in the literature, since 
it provides immediate information on the strength of heating 
by conduction in different clusters. Starting from Eqs. |3 - 16] 
one can estimate whether heat conduction alone can be efli- 
cient in quenching the cooling flow by setting both H and M 
equal to zero. For a sample that includes all the objects stud- 
ied in the present work, K04 have explored this possibility and 
concluded, in agreement with other studies (Voigt et al. 2002 
Voigt & Fabian E004L that heat conduction alone is insufficient 
to balance radiative losses in CF clusters (with the exception of 
MKW 9 and A 399, see also discussion below), in particular 
in their central regions. These results and the growing obser- 
vational evidence of AGN-ICM interaction at the center of CF 
clusters motivate including in the model an additional heating 
term provided by a central AGN. 

4.2. Conduction and AGN heating: zero mass-dropout 

Motivated by the results of RB02, who find that the efferves- 
cent heating model does not require any mass dropout, we first 
present results for M = 0. As this case is the most justified, we 
present its predictions in more detail. In addition, as we show 
below, small mass-inflow/outflow rates do not modify the re- 
sults substantially. For M = 0, we vary fc between and 1/3 
and fit the extra heating curve H (Eq.|6j using the heating func- 
tion i/'^CAf (^£q i\ ]sjo constraints are imposed on the AGN 
luminosity L, but the best-fit values for ro are searched in the 
interval from to the cooling radius RcooU which is computed 
for a cooling time of 15 Gyrs and is listed in Table|3] 

In Fig. [Owe illustrate results for the cluster A 1795 and 
discuss them in the following, since the model outcome for this 
object highlights important features also found for most of the 
clusters in the sample. Since no mass dropout is present, e* is 
zero in Eq.|6land the extra heating H - b - Econd only depends 
on the conduction efficiency f^. Since = 1 /3 is the maximum 
value we consider, it corresponds to the maximum energy yield 
by heat conduction from the outer parts of the cluster. From a 
visual inspection of Fig. [2 it is clear that heat conduction is 
not able to lower the extra heating in the outermost bins. As a 
consequence, if the extra heating curve is fitted with Eq.[^ the 
resulting best-fit parameter ro (the scale radius where the bub- 
bles start rising in the ICM) is unphysically large (ro = 176 kpc 



> ^cooi, see Table |5} and the time-averaged AGN luminosity 
L = 3.3 X 10"^^ erg s"' is also quite large. The latter solution is 
not taken into account since, as stated above, the best-fit value 
for ro is searched in the interval from to R^ooh but is used here 
as an example. 

On the other hand, if the effervescent heating is assumed 
to be efficient in the region only within the cooling radius, 
one can notice that heat conduction is effective in reducing 
the extra heating so that the resulting extra heating curve is 
monotonically falling with radius. In this case one therefore 
expects that the extra heating supplied by the raising relativis- 
tic bubbles must be distributed over smaller distances and that 
the total AGN energy output is lower. This is indeed reflected 
in the much different best-fit parameters ro = 14 kpc and 
L = 8 X 10^"^ erg s"' . This feature indicates that, as expected, the 
effervescent heating model strongly depends on how much and 
at which radial distance heat conduction is efficient. Because 
of this, we performed, for each cluster and each analytic func- 
tion used to model the gas density and temperature, fits using 
the whole observed radial range (i.e., < r < /?T.max) and the 
radial range delimited by the cooling radius. 

Another common feature is the effect of the variation of the 
conduction efficiency. If /c = the extra heating curve is sim- 
ply equal to the gas emissivity and the increase in /c from to 
1/3 gives a decrease in the extra heating curve from the emis- 
sivity curve to the data points marked in Fig. ^ The increase 
in /c should hence lead to a decrease in the AGN energy re- 
quirement. In fact we find that, within the cooling radius, both 
time-averaged AGN luminosity and inner cutoff radius ro de- 
crease monotonically with increasing f^. This opposite effect is 
seen if the fits are performed over the whole radial range: an 
increase in /c leads to an increase in both L and ro, again show- 
ing the inadequacy of applying the model over the whole radial 
range. 

For the majority of the objects in our sample, we find that 
the effervescent heating model provides results when the radial 
range used to fit the extra heating curve with Eq.[nis not sim- 
ply the whole observed range but the radial range inside the 
cooHng region. In fact for NGC 533, Virgo, A 262, MKW 9, 
A 4059, and A 496, only fits performed over the radial range 
inside the cooling region converge. The reason for this lack of 
convergence of the fits is, as in the case of A 1795 illustrated 
in Fig. [2 that heat conduction is not able to reduce the extra 
heating in the outermost bins. On the other hand, a reduction of 
the extra heating is present in bins close to the cooling radius. 

This is a welcome feature for the effervescent heating 
model, because it implies that heating from the central source 
is not only more viable in terms of total energy output, but also 
that the energy can be distributed from the center outwards up 
to relatively small distances. For Sersic 159-3, 2A 0335-1-096, 
Hydra A, and A 1835 heat conduction is low even for the max- 
imal conduction efficiency fc - 1/3. Therefore a characteris- 
tic drop of the extra heating close to the cooling radius is not 
present in these objects and hence no difference is seen between 
the best fits performed using different radial ranges. The low- 
ering of the extra heating close to the cooling radius is also not 
seen in MKW 3s, A 2052, and A 3112, because heat conduc- 
tion is quite efficient only at the cluster center. While in A 3 1 12 
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Fig. 1. Energy requirements in A 1795. The plasma emissiv- 
ity is shown, for simplicity, by the dashed straight lines join- 
ing the emissivity values computed in bins and heating due to 
thermal conduction by the dot-dashed line. The mass flow rate 
is set to zero in this model, so the extra heating curve (filled 
squares) is determined by f^x Spitzer. Here /c = 1 /3, the max- 
imum value allowed in the models. The three bins inside the 
cooling radius (vertical line) alone are best fitted by an inner 
cutoff radius ro = 14 kpc and time-averaged AGN luminos- 
ity L = 8.0 X lO'^'^ergs"^ For the same model (M - and 
/c - 1/3), all six bins are best fitted by an unphysical ro = 179 
kpc and very high AGN luminosity L = 3.3 x 10"*^ erg s"' . The 
two solid curves are the best-fit functions for the two cases. See 
text for details. 

heat conduction is low and therefore the best-fit parameters for 
the two radial ranges are almost identical, for MKW 3s and 
A 2052 we obtain slightly higher values for both ro and L when 
the whole radial range is used instead of the cooling region. 

For these reasons, only results from fits to the extra heating 
curves within the cooling radius will be discussed. For A 399, 
no results are obtained for either of the two radial ranges. The 
reason for the lack of convergence of the fits is that in A 399 
heat conduction is very efficient, owing to its relatively high 
temperature and very steep temperature profile (see also K04, 
Table 1 1). Conduction can balance radiative losses completely 
for fc as low as 0.03. For lower values of the conduction effi- 
ciency the best-fit values for ro are always larger than the cool- 
ing radius Rcooi = 43 kpc and therefore not taken into account. 

The results for the zero mass dropout models are summa- 
rized in Table|3] For each value of fc the 1 cr errors on the best- 
fit parameters are very small compared to the range spanned by 
the best-fit parameters for the different models and therefore 
are not reported. In addition, it is important to notice that when 
fits are restricted to the cooling region, the number of bins we 
are using in the fits is quite low (3 for 9 objects, 4 for 4 ob- 
jects, and 5 for 3 objects), and although it is sufficient for our 
scope, it should be noticed that for more complex models (e.g., 
models with yb as an additional free parameter) the number of 
bins within the cooling radius should be increased if possible. 
While for most of the clusters we find a solution for L and ro 
for every value of < fc < 1/3, for some objects the fits con- 
verge only for models with f. in a narrower interval (/c™"'-/c"^'', 
see Table|3j. In particular, in clusters where conduction is high 



(MKW 9 and Perseus), solutions are found for conduction ef- 
ficiencies substantially lower than the maximum allowed value 
fc = 1/3. In Sersic 159-3, 2A 0335+096, Hydra A, A 3112, 
and A 1835, the conductivity yield is low and, as expected, dif- 
ferent models (i.e., different values of fc) give almost identical 
results. 

It is crucial to highlight the trend of the best-fit parameters 
of the different models with model parameter f-. For most of 
the clusters, the model with /c = /f"^" is the one where val- 
ues for both L and ro are lowest and for /c - /™" the high- 
est (L"in.™'< and r™"'""" in Table|3. In addition L and ro vary 
monotonically with fc within these Umits. The clusters that ex- 
hibit this trend are labeled by an asterisk in column trend of 
TableE For MKW 3s, A 2052, and A 31 12 the trend just de- 
scribed is reversed. In these objects heat conduction lowers the 
extra heating curve especially in the cluster center This implies 
more and more flattening of the extra heating curve at the cen- 
ter with increasing fc, which is finally reflected in the increase 
in both ro and L. In Table|3lwe also list the best-fit parameters 
L and ro for the best model (fc), which is chosen by simply 
comparing the;^^^ of the best fit for the different models. 

4.3. Conduction and AGN heating: models with mass 
dropout or outflow 

Although RB02's results indicate that an equilibrium state with 
no mass dropout is usually achieved, the effervescent heating 
model incorporates the possibility of gas inflow or outflow. It is 
therefore interesting to see how much our results are changed 
by including gas inflow/outflow (negative/positive mass flow 
rates), since this implies a positive/negative work (per unit vol- 
ume) done by the system (e* in Eq.|6}. It is expected that neg- 
ative mass-flow rates reduce the energy required by the central 
source to quench cooling, but the question is how large this 
effect is. 

We apply the effervescent heating model as done in 
Sect. 14.21 but we go on to explore models with -M™™/10 < 
M < -i-M™'"'/10. When possible, the maximum mass flow rate 
^max ggj. gq^^J jjjg maximum of the mass deposition rates 
given in Table 5 in Peterson et al. ( I2003t . Otherwise M™'"' is 
set equal to the maximum value that, at /c = 0, still gives 
positive values for H inside the cooling radius. The values of 
the maximum mass flow rates M™^" are listed in Table |3 This 
choice of the upper and lower bounds for M is dictated by the 
fact that we find the effect of a fixed mass inflow/outflow rate 
for all the clusters to be drastically different from cluster to 
cluster In addition, taking RB02's results into account, low 
mass flow rates should be present when an equilibrium state 
(which is assumed here) is achieved. The analysis of models 
with -M^^/IO < M < +M'™''/10 is therefore a good com- 
promise. We do not present the results from fitting the mod- 
els extensively. More important we find that, as expected, the 
effect introduced by mass inflow or outflow is to reduce the 
time-averaged luminosity L for inflows and enhancing it for 
outflows. From the models we have studied and, of course, for 
those where the fits converge in the L - ro plane, we report min- 
imum and maximum values for the time-averaged luminosity 
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Table 3. Cooling radius 7?cooi and the results for the effervescent heating model with no mass dropout (M = 0). For each model, 
i.e. for a fixed conduction efficiency < /c < 1/3, the extra heating curve H derived from Eqs.|3-|6lis fitted using Eq. 
over the radial range < r < R^oh and the best-fit values for L and ro are derived. Only models with between /™'" and Z™" 
can be fitted. The corresponding minimum and maximum values of the time-averaged luminosity L and the inner cutoff radius 
ro are given. For most of the clusters, both best-fit parameters for time-averaged luminosity and inner cutoff radius decrease 
monotonically with increasing /c. These are labeled by an asterisk in the column trend. L and ro are the f/me-averagec/ luminosity 
and the inner cutoff radius for the model with the smallest reduced and fc is the corresponding conduction efficiency. The 
minimum and maximum values of the time-averaged luminosity z^™"'™" for finite integration limits (see Sect. 14. 5> are also 
listed. No result is found for A 399. 
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Table 4. Results for the effervescent heating model with mass dropout or outflow. The minimum and maximum values of the 
time-averaged luminosity L and the inner cutoff radius ro are computed from the results of models with conduction efficiency 
< /c < 1/3 and -M^^/IO < M < +M™''/10. The minimum and maximum values of the f/me-averagec/ luminosity z^™"-™" 
for finite integration limits (see Sect. 14. 5> are also listed. As for the M - Q models, no result is found for A 399. 
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L and the inner cutoff radius ro in Table |3 A comparison with 
the values for the model with no mass dropout (Table|3} shows 
that the broadening of the range for the best-fit parameters is 
not large. 



4.4. Effects of the temperature and density modelling 

In order to investigate the effects introduced by the choice of 
fitting functions for the density and temperature profiles, we 
performed the same analysis presented up to here (i.e., using 
Eq. n and Eq. |3} but using the remaining possible combina- 
tions of parametrizations given in Sect. 12.11 We find that for 
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a fixed model (i.e., for fixed M and fc) the use of different 
parametrizations introduces small differences in the estimated 
best-fit parameters. These are smallest for the models with- 
out mass dropout or outflow and in particular for those clus- 
ters in which heat conduction is inefficient. For Sersic 159-3, 
2A 0335-H096, Hydra A, and A 1835, the difference in both L 
and ro is less than 3% and for the rest of the clusters less than 
10%. The differences are largest when the double yS-model and 
Eq. 121 are used to parametrize gas density and temperature for 
models with M + Q. Also in this case, the difference is smallest 
for the clusters in which heat conduction is not efficient (typ- 
ically less then 7% for both best-fit parameters L and ro). For 
the rest of the clusters, the difference can be as high as 18%, but 
generally less than 14%. Taking other parametrizations of the 
temperature profile like those used in Churazov et al. (2003 1 
and Dennis & Chandran J2005t leads to a change of at most 
10% for both best-fit parameters L and ro for a fixed model. 
Most important, the trends between the AGN parameters L and 
ro and the model parameters fc and M described in Sects. l4!2l 
and l4.3l remain unchanged for every combination of gas density 
and temperature parametrizations used. We therefore conclude 
that the values given in Tables |3 and |3 give quite robust esti- 
mates for the time-averaged luminosity L and the inner cutoff 
radius ro. 



4.5. Effects of the integration limits 

Since most of the models in which radiative cooling is bal- 
anced by some heating source fail due to insufficient energy 
yield, it is important to investigate to what extent our values 
for the time-averaged luminosity of the central AGN are over- 
estimated. When fitting the extra heating function (Eq. II 1> to 
the extra heating curve, one directly derives the best-fit val- 
ues for ro and the normalization of the extra heating function, 
i.e. LjAnq. The time-averaged luminosity L is then estimated 
using Eq. ^] The values of L reported so far were computed 
by setting the integration limits as given in Eq. but as al- 
ready noted in Sect.|3j the upper integration limit should be re- 
placed. As reported in Sect. 14.21 the effervescent heating model 
provides results for almost all the clusters when applied to the 
cooling region. As a consequence we set the upper integration 
limit in in Eq.^jequal to the cooling radius Rcoo\- 

The effect of this replacement is to reduce the time- 
averaged luminosity of the central AGN. The effect of this 
simple but effective correction is shown in Tabled where we 
list the minimum and maximum time-averaged luminosities 
L™" and L™" computed using this finite integration limit for 
the models with no mass dropout. The same quantities for the 
model with mass dropout are listed in Table|3 Comparing these 
values with the ones computed in the previous sections shows 
that the reduction in the time-averaged luminosity can be al- 
most one order of magnitude, hence a finite integration limit in 
Eq.[Olshould be used when estimating the effervescent heating 
best-fit parameters. 



5. Discussion and conclusions 

We have used deprojected radial density and temperature pro- 
files of a sample of 16 nearby CF clusters observed with XMM- 
Newton (Kaastra et al. 2004 1 to test whether the effervescent 
heating model (Ruszkowski & Begelman l2002l l can satisfacto- 
rily explain the structure of CF clusters. The effervescent heat- 
ing model incorporates both heat conduction and AGN feed- 
back as heating sources, which is of great interest because of 
the complementary nature of these two processes. For each 
cluster, we derived the required extra heating as a function 
of cluster-centric distance for various values of the unknown 
parameters M (mass deposition rate) and fc (conduction ef- 
ficiency). We fitted the extra heating curve using the AGN- 
heating function proposed by Ruszkowski & Begelman J2002> 
and derived the AGN parameters L (the time-averaged lumi- 
nosity) and ro (the scale radius where the bubbles start rising in 
the ICM). 

For models without gas mass dropout or outflow (M - 0), 
we find: 

- for only one object (A 399) we do not find any solution for 
the effervescent heating model because heat conduction is 
very efficient; 

- for 4 clusters (Sersic 159-3, 2A 0335+096, Hydra A, and 
A 1835) the conductivity yield is extremely low and, as ex- 
pected, different models (i.e., different values of fc) give 
almost identical results; 

- for 3 objects (MKW 3s, A 2052, and A 3112), we find that 
heat conduction plays an important role only at the cluster 
center and that, as a consequence, the trend between the 
fitted AGN parameters and conduction efficiency is not the 
one expected if conduction and AGN heating are assumed 
to cooperate; 

- for the remaining 8 clusters (NGC 533, Virgo, A 262, 
MKW 9, A 4059, A 496, A 1795, and Perseus), conduc- 
tion and AGN heating are found to be cooperating. 

We have studied models with mild gas mass dropout of outflow 
(-M™V10 < M < +M™V10, with M™"" being the mass de- 
position rate required by the standard CF model) and find that: 

- the conclusions listed above for the case of no mass 
dropout/outflow are also valid in this case; 

- as expected, the effect introduced by mass inflow or outflow 
is to reduce the required AGN heating for inflows and to 
enhance it for outflows; 

- this implies that the ranges of allowed AGN parameters L 
and ro is broader than in the case M = 0, but we show that 
the broadening is not large. 

We find that our results are not sensitive to the choice of fitting 
functions used to model gas density and temperature profiles 
and that the time-averaged AGN luminosities required to bal- 
ance radiative losses are substantially reduced if the fact that 
the AGN deposits energy within a finite volume is taken into 
account. 

Since we do not find any solution with the effervescent 
heating model for only one object (A 399), we conclude that 
the model provides a satisfactory explanation of the observed 



10 



Piffaretti & Kaastra: Effervescent heating in CF clusters 



structure of CF clusters. As pointed out by Begelman (I2004> . 
AGN-heating is dominant in the final stable state of the clus- 
ter modelled by Ruszkowski & Begelman (2002i. Conductive 
heating dominates only at the beginning of the evolution and 
might be essential only for the stability of the model. Therefore, 
models with high heat conduction are not preferred. 

On the other hand, since there is evidence that AGN heat- 
ing alone is not able to quench CFs (Brighenti & Mathews 
120021 Zakamska & Naravan I2003> . it is fair to assume that 
thermal conduction, although operating at different rates from 
cluster to cluster, must play an important role as a heating 
mechanism. While this is found for 8 clusters in our sample 
(NGC 533, Virgo, A 262, MKW 9, A 4059, A 496, A 1795, and 
Perseus), we have shown that conductive heating is either com- 
pletely unimportant in 4 clusters (Sersic 159-3, 2A 0335+096, 
Hydra A, and A 1835), too high for one object (A 399), or high 
enough to play an important role but peaked at the cluster cen- 
ter in 3 clusters (MKW 3s, A 2052, and A 3112). Therefore, 
if we assume that AGN and conduction heating must be co- 
operating effectively, the model does not provide a satisfactory 
explanation for half of the objects in the sample. 

These findings prompt us to posit that, at least for these ob- 
jects, the description of their thermal structure through a steady 
state solution of the thermodynamic equations is not viable and 
that we are observing them at an evolutionary stage far from 
equilibrium. A clearer picture can, of course, be achieved by 
studying a much larger sample using the procedure employed 
in this work. 

For the clusters that are well described by the effervescent 
heating model, the derived best fit parameters ro and L can be 
compared with results obtained from complementary observa- 
tions. In the effervescent heating scenario, it is supposed that 
bubbles are deposited by the central radio source. Although the 
nature of this process is not well understood, an interesting pos- 
sibility is that bubbles are generated through the interaction of 
the radio jets with the surrounding ICM. In the framework of 
the effervescent heating model, the AGN parameter ro fixes the 
scale radius where the bubbles start rising in the cluster atmo- 
sphere and heat the ICM. One then expects the AGN parameter 
ro to be larger than the jet extension and of the same order of 
magnitude. Since the observed jet extension rjet is at most as 
large as the true jet size due to projection effects, the same re- 
lationship is expected between ro and rjet. 

The size of the AGN jet has only been measured for 10 
of the 15 clusters that are well described by the effervescent 
heating model. For these 10 objects, we show the comparison 
between ro and rjet in Fig.|2 Even though the range of allowed 
values for ro is quite large, the comparison shows that, for most 
of the clusters, the extension of the AGN jet and the scale radius 
ro are within the same order of magnitude and that the latter is 
in general larger than the jet extension. This implies that the 
possibility of bubbles being generated through the interaction 
of the radio jets with the ICM cannot be excluded. 

A comparison between the derived luminosities L with ob- 
served luminosities is by far more difficult. In fact, in the frame- 
work of the effervescent heating model, the derived AGN lu- 
minosity is a time-averaged total AGN power. Since the ra- 
dio power is a poor tracer of the total power (Eilek .2004j a 
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Fig. 2. Comparison between the observed extension of the 
AGN jet rjet and the best-fit AGN parameter ro. For ro we 
plot the whole range of allowed values for the models with no 
mass dropout (see Table O, and the extension of the AGN jet 
has been taken from the literature: Virgo (Young et al. I2002> . 
A 262 (Blanton et al. 2004i, 2A 0335+096 (Sarazin et al. 
[T995t , A 2052 (Blanton et al. HM^, A 4059 (Taylor et al. 
[T994l l. Hydra A (McNamara et al. l2000t . A 496 (Markovic et 
al.l2004l l. A 31 12 (Takizawa et al. 2003 1, A 1795 (Ettori et al. 
I2002t . and Perseus (Fabian et al. .2002j . The solid line indicates 
the locus of equality. 

fair comparison is possible only for objects with estimated to- 
tal jet power. At present, estimating the total jet power is un- 
fortunately possible for only one source: M 87 in the Virgo 
cluster (Owen et al. I2000t . In this case, our estimate of L 
(see Table O agrees well with the inferred total jet power, 
~ 3 - 4 X 10"*^ ergs"^. A similar conclusion is presented in 
Ghizzardi et al. J^nO^t . 

Even though we used the largest sample of CF clusters with 
well-measured deprojected gas density and temperature pro- 
files available at present, it is clear that its size is not large 
enough to draw statistically significant conclusions on the via- 
bility of the effervescent heating model. Therefore, considering 
the simplicity of the procedure, any future work with the aim 
of extending the analysis to a much larger sample will provide 
vital information on this issue. 
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